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Abstract 

We review the calculation of the prompt lepton flux, produced in the atmosphere 
by the semileptonic decay of charmed particles. We describe side by side the inter- 
mediary ingredients used by different authors, which include not only the charm 
production model, but also other atmospheric particle showering parameters. After 
evaluating separately the relevance of each single ingredient, we analyze the effect of 
different combinations over the final result. We highlight the impact of the prompt 
lepton flux calculation upon high-energy neutrino telescopes. 
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1 Introduction 

Very-high-energy (above 1 TeV) neutrino astronomy is currently a subject of 
great interest, promising to expand our observational range of the Universe in 
an unique way[l]. Such energetic neutrinos may carry information from the 
sources of the highest energy phenomena ever observed in cosmic rays, pos- 
sibly coming from active galactic nuclei (AGN) or gamma ray bursts (GRB). 
They may probe the early stages of the Universe and its farthest distances. In 
addition, they will contribute to the search of weakly interacting massive par- 
ticles (WIMP), supernova explosions, monopoles, besides the discovery poten- 
tial for new physics, which we can even not imagine yet. Neutrino telescopes 
under development, like the Antarctic Muon and Neutrino Detector Array 
(AMANDA) and the experiment at Lake Baikal, are already operational, and 
producing their first results[2,3]. In addition, great activity is planned for the 
near future [4]. 

Aside from these perspectives, the operation of a neutrino detector at energies 
above 1 TeV poses challenging difficulties. One of the major limitations in 
the detection of a cosmic high-energy neutrino (from galactic or extra-galactic 
origin) is the background from atmospheric muons and neutrinos (produced 
by the interaction of high-energy cosmic rays in the atmosphere) . 

The source of the atmospheric neutrino background changes with energy, in 
a way governed by the critical energy e cri u c of the parent particle. This is 
the energy for which the decay and interaction lengths are equal. Above this 
energy the parent particle is likely to interact or be slowed down before decay- 
ing into a neutrino. The critical energy is calculated in terms of the particle 
rest energy mc 2 , the mean life r and, by adopting the isothermal atmosphere 
approximation, a scale constant h Q [5]: 

mc 2 

^critic i^o 
CT 

Table 1 summarizes several particle properties, derived from the Review of 
Particle Physics [6] . Comparing the critical energies we observe that muon de- 
cays contribute substantially to the atmospheric lepton flux only up to a few 
GeV's, while the decays in flight of charged pions and kaons are still impor- 
tant up to about 1-10 TeV. They give rise to the "conventional" atmospheric 
lepton flux. Above this energy range, the semileptonic decay of very-short 
lived charmed particles (mainly D-mesons and A+-hyperons) is the dominant 
source, despite their low production rate. The main contribution comes from 
the decay modes 

D — > K + ji + v and A c — > A c + /x + v. 
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Muons and neutrinos thus generated are called "prompt leptons", and they 
exhibit a flatter (and thus harder) energy spectrum. The lack of precise in- 
formation on high-energy charm production in hadron-nucleus collisions leads 
to a great uncertainty in the estimate of the leptonic flux above 100 TeV. In 
addition, different authors do not use the same atmospheric particle shower- 
ing routines, turning the comparison even more difficult. The predictions of 
resulting fluxes span up to four orders of magnitude! 

It is our purpose to bring some light to the forum of prompt lepton fluxes, 
describing side by side the many ingredients of the calculation. After evaluat- 
ing separately the relevance of each single shower parameter, we analyze the 
effect of different combinations over the final result. We investigate the fluxes 
of \i, and v e , leaving the case of r and v T for a further work. 

Table 1 

Particle data summary. 



Particle 


Elementary 


mc 2 


cr 


£■ critic" ^ 


Branching 




contents 


(MeV) 




(GeV) 


ratio £j( 2 ) 


D+,D- 


cd, cd 


1870 


317 um. 


3.8 x 10 7 


17.2 % 


D°,D° 


cu, cu 


1865 


124 jum 


9.6 x 10 7 


6.8 % 


Df,D~ 


cs, cs 


1969 


149 ^m 


8.5 x 10 7 


5.2 % 


A+ 


udc 


2285 


62 um 


2.4 x 10 8 


4.5 % 


u + ,u~ 


lepton 


106 


659 m 


1.0 


100 % 


7T + , 7T~ 


ud, ud 


140 


7.8 m 


115 


100 % 


K + ,K- 


us, us 


494 


3.7 m 


855 


63.5 % 


A° 


uds 


1116 


7.9 cm 


9.0 x 10 4 


0.1 % 



(1) According to Eq. (1), with h = 6.4 km. 

(2) For inclusive decays yielding leptons. 



2 Calculation of the prompt lepton flux 

The calculation of the prompt lepton flux has been carried out in the past 
(see, e.g. Ref.[7] and references therein), mainly with the purpose to investi- 
gate the effects of choosing a different charm production model. The outline 
of the calculation is basically the same in all of these works. We start from 
the primary cosmic ray flux at the top of the atmosphere, with a composition 
supposed to be dominated by protons, and evaluate the flux of nucleons at any 
atmospheric depth. Those nucleons interact with the nuclei of air to produce 
secondary particle showers. For energies above a few TeV, the secondary par- 
ticles of interest to be followed are the charmed hadrons, for they will be the 
main source of high energy atmospheric leptons. We will, therefore, consider 
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the contribution of mesons D ± ,D° and D° (comprising in the same notation 
the overall contribution of D and D* mesons), of the strange Df and of the 
A+-hyperon. 

Because the production rate of the mesons D s is relatively lower (about 20% 
of the D production cross section), some authors neglect their contribution, 
although their critical energy and branching ratio are comparable to those of 
the other charmed particles (see Table 1). When calculating the flux of r and 
u T , which we do not consider in this analysis, the role of D s turns out to be 
mostly important [8]. 

It is straightforward to calculate the flux of charm particles at any depth, 
and they will promptly decay yielding electrons, muons and neutrinos. We 
may integrate the flux for all possible charmed parent production and decay 
depths, and for all possible production and decay energies, leading to the flux 
of a chosen lepton at a given depth and energy. 

For detailed calculations we refer to Refs. [7,9,10], and we follow the notation 
of the latter, to present here only main results. Let's write the primary cosmic- 
ray spectrum as a power law in energy: 

$ N (E N ,x = 0) = N E N h+1 \ (2) 



where &n{En, x), given in (GeV.cm 2 .s.sr)~ 1 , is the differential flux of nucleons 
with energy En, in GeV, and x is the slant depth penetrated by the cascade, 
measured in g/cm 2 from the top of the atmosphere (x = 0) downward along 
the direction of the incident nucleon. The constant N Q is the amplitude, or 
differential spectrum normalization; and 7 is the spectral index, or slope of 
the integral primary spectrum. 

After developing to a certain depth x, the nucleonic flux is given in terms of 
A N , the nucleonic attenuation length [5]: 

$ N (E N , x, 6) = N Q E N { " +1) e^ A -. (3) 



The resulting flux of secondary particles of type-i (i = D ± , D°, D°, Df, A+) is 
calculated by convolution of the nucleonic flux with the production spectrum 
of secondary particles: 

MEi,x,0) = KiiEu-y) Jdx> (^J expj-^-^ - £-| , (4) 
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defining 

^critic 



Ei cosO' 



where the dependence in the zenith angle holds for 9 < 60°. For higher zenith 
angles the curvature of Earth must be taken into account. Both the nucleonic 
attenuation length A^r and the charmed particle interaction length Aj are given 
in units g/cm 2 . The production spectrum of charmed particles, weighted by 
the primary nucleonic spectrum, is written as: 



N 

K l (E l , 1 )= I dE N E 



A 



N 



J N 



-(7+1) 



i) dW Ni (Ei,E_ 



dEi 



= ^Ei^Z m { 1 ). 
an 



(5) 
(6) 



In this notation, dW Nl /dEi denotes the energy distribution of secondary par- 
ticles, and represents the probability that a particle of type-i with energy Ei 
is created in the interaction of an incident nucleon N of energy E^ with an air 
nucleus. It is directly related to the inclusive cross section for secondary parti- 
cle production. Eq. (6) is obtained assuming a mild energy dependence for the 
nucleonic interaction length \ N , and defining Z^^), the particle production 
spectrum- weighted moment [5], also called production "Z-moment": 




where xf is the Feynman variable xf = Pl/Pl™, with pi as the produced 
particle longitudinal momentum. At the high-energy limit the Feynman-rc also 
represents the ratio of the final particle energy to the incident particle energy, 
xp = Ei/E N (beware confusion with atmospheric depth x). 

In order to evaluate the flux &i(Ei,x, 9) of leptons {I — n or u), with energy E\ 
and zenith angle 9 at depth x, we need to fold the energy distribution df l /dEi 
of the produced lepton with the spectrum Di(Ei,x" ,9) of decaying parents, 
for any decay depth x" and any available parent energy Ei\ 

jpmin 

X i 

$i(E h x,6) = jdx" J dEi -jL Di{E h x",9), (8) 

pmax ' 



with 
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Di(Ei, x", 6) = Bij ®i(E h x", 6), (9) 

x" x"cos0 Ei 
di = — = (10) 

f] ^critic 



where Bi is the branching ratio yielding leptons in the parent-i decay (see 
Table 1), and di is the particle- i decay length. 

The muon and neutrino production energy distributions df 1 / dE\ used in Eq. (8) 
are given by the semileptonic three-body decay phase space integrals, ob- 
tained from kinematics considerations^, 11]. Some authors[12-14] define a de- 
cay Z-moment in analogy to the production Z-moment, Eq. 7, just replacing 
dW Nl / dEi by df 1 jdE\. Making use of this definition, it is possible to write an 
approximate solution (valid for energies E\ < e C riUc) exploring the fact that the 
critical energy for charmed particles is very high. Nevertheless, in the present 
work we will use the complete solution for calculating the prompt lepton flux, 
given by the set of Equations (4) to (10). 



3 Showering Process 



Once the calculation is established, the next step is to choose the ingredients 
that characterize the showering process in the atmosphere. We will browse 
through the literature to extract different parametrizations to be compared. 
The main parameters to define are: the primary spectrum normalization and 
slope, the nucleonic and charm interaction lengths, the nucleonic attenuation 
length and the charm production Z-moment. 



3.1 Primary spectrum 



The primary cosmic ray flux at the top of the atmosphere, Eq. (2), can be 
rewritten as to incorporate the effect of the change in slope ( "knee" ) observed 
in the energy spectrum, at energy E knee : 

<£> N (E N ,x = 0)=m E N ^ 1+1 \ E N < E knee (in GeV); 

= N 2 E~ N {12+1 \ E N > E knee (in GeV). (11) 

Table 2 indicates some typical values for the parameters in Eq. (11), used 
by different authors. Lipari[12] quotes a parametrization consistent with both 
the JACEE balloon borne experiments[15] and the values given by Gaisser[5], 
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Table 2 

Frimary cosmic-ray fluxW 



Label 


E <C Ek nee 


Eknee 


E > Eknee 




Ni 


71 


(GeV) 


N 2 72 


Lipari 


1.70 


1.70 






Akeno 


1.35 


1.62 


4.67 x 10 6 


630 2.02 


Bugaev (F) 


1.02 


1.62 


1.9 x 10 6 


323 2.02 


Bugaev (D) 


1.02 


1.62 


5.2 x 10 5 


193 2.02 


TIG 


1.70 


1.70 


5 x 10 6 


174 2.00 


(1) According 


; to Eq. 


(11), 







with Ni and N 2 given in units (GeV.cm 2 .s.sr) 1 . 

adopting a single slope, because his analysis is mainly aimed at energies below 
the knee. The AKENO experiment [16] obtained a description of the primary 
spectrum, covering the knee region, from data on size spectra of electrons 
and muons at high energy. Bugaev et al. [17] use a semiempirical model which 
takes into account detailed chemical composition of the primary spectrum, 
translated here in terms of Eq. (11). They propose two options (Model F and 
Model D), differing on the hypothesis responsible for the change in slope at 
the knee. Thunman, Ingelman and Gondolo (TIG) [13] follow the JACEE 
trend below the knee. Above the knee they adopt the same slope as Volkova 
et al.[9], for which no normalization constant is reported, due to their interest 
limited to flux ratios. Figure 1 displays the energy spectra corresponding to 
those parametrizations. 



3.2 Interaction and attenuation lengths 



The nucleonic interaction length, A at, represents the mean free path of nucleons 
in the atmosphere (given in g/cm 2 ). It is related to o-^~ atT ', the total inelastic 
cross section for collisions of nucleons (N) with air nuclei, through the relation 

A.m p 24100 



where we used the average atomic number for air nuclei A—1A.5, m p is the 
proton mass and o^~ mr must be given in mb. There are different parametriza- 
tions to the inelastic A^— air cross section as a function of energy. Some au- 
thors considered it to be constant [9, 12]; others to be rising with energy: as a 
power-law[16], as a logarithmic dependence [7, 17], or as a log-squared depen- 
dence [18, 19]. 
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Fig. 1. Comparison of primary cosmic-ray energy spectra, as given by different 
parametrizations (see Table 2). The primary flux is multiplied by E N , so that the 
structure at the knee can be better appreciated, and the different parametrizations 
compared. 

The attenuation length, An, which governs the exponential decay of the nu- 
cleonic flux with increasing depth, see Eq. (3), represents the net effect in 
the interplay between interaction losses and regeneration of the number of 
nucleons in the cascade development, and is given by 



A 



A 



N 



N 



(i - z NN (ij) 



(13) 



The interaction length Xn dictates the losses and the nucleon-to-nucleon spec- 
trum-weighted Z-moment Znn{i) accounts for the survival rate of nucleons. 
Z NN is calculated in analogy to Eq. (7), with the outgoing particles be- 
ing the regenerated nucleons. Some authors use an approximately constant 
value[7,12,19], others adopt Feynman- scaling with a change in value at the 
knee energy [5, 18], while in Ref. [13] it is assumed the violation of Feynman- 
scaling. 

The charm particle of type-i may have an interaction length Aj in the atmo- 
sphere calculated analogously as for the nucleons, Eq. (12), substituting a^~ air 
by cr*~ aM '. As in the parametrization of the nucleonic inelastic cross section, we 
find authors adopting a charm cross section which is constant [20], or which 
increases with energy either as a power-law [21], or logarithmically [7]. 
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3.3 Charm production spectrum-weighted moments 



The key information for the evaluation of the prompt lepton flux is the be- 
havior of the charm spectrum-weighted moments, given by a specific charm 
production model. Among these models, we select three different ones, to be 
compared in the present study. 

QGSM: Quark Gluon String Model, a semiempirical model of charm pro- 
duction based on the non-perturbative QCD calculation by Kaidalov and 
Piskunova[22], normalized to accelerator data, and applied to the prompt 
muon calculation by Volkova et al.[9}. 

The Z-moments are calculated numerically from Eq. (7), with 

dW m _ cf NA dU 



dx F a^ n air dxp 
where the particle production spectrum is parametrized by: 

d fi L08 r , 

(1 — Xp) for D production, 



dxp Xp 

= 1.4 (1 - x F f A for A c production. (15) 

For the total inelastic cross section, of n ~ air , we used a parametrization with a 
log 2 energy dependence [18]. Following Ref.[9], a mild log energy dependence 
is assumed for the inclusive cross section of charm production, cf NA , and the 
production of D s is neglected. Figure 2 displays the curves of Z-moments 
for the different charm components. 

RQPM: Recombination Quark Parton Model, a phenomenological non-per- 
turbative approach, taking into account the contribution of the intrinsic 
charm to the production process, in which a cc pair is coupled to more than 
one constituent of the projectile hadron, as described by Bugaev et al. [7, 17]. 

Supposing Feynman-scaling to hold (RQPM-FS), the charm production 
Z-moments are given simply by: 

Zm{l) = Z y = constant, (16) 

with parameters for different particles shown in Table 3. Assuming the scal- 
ing violation (RQPM-SV), the parametrization turns out to be: 




Zm{l) = Z, -f , (17) 



where £ = 0.177 — 0.05 7. The parameters are also given in Table 3. We 
obtained the all-charged D-meson Z-moment (i = D^) by averaging the 
individual conjugate moments, the same being done for neutral D-mesons, 
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Fig. 2. Charm Z-moment components for QGSM, as a function of energy. 
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Fig. 3. Charm Z-moment components for RQPM, considering respectively Feyn- 
man-scaling (FS) and scaling violation (SV), as a function of the energy. 

{% = D,D°). Ref. [7] does not take into account the contribution of D s 
mesons. A comparison of the resulting Z-moments, with and without scaling, 
is provided in Figure 3. 

pQCD: Perturbative QCD, as calculated by TIG [13], using the Monte Carlo 
program PYTHIA[23], explicitly evaluating the charm production, up to 
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Table 3 

RQPM Z-moment parameters Z 7 



Label 


7 


e 


E 1 (GeV) 


Z 7 (D±) 


Z 7 (D°,D°) 


Z 7 (A+) 


RQPM-FS 


1.62 






4.88 xl0~ 4 


4.73 xlO~ 4 


4.36 xlO" 4 




2.02 






3.14 xl0~ 4 


3.09 xlO" 4 


2.95 xlO" 4 


RQPM-SV 


1.62 


0.096 


10 3 


5.55 xlO" 4 


5.35 xlO" 4 


4.9 xlO" 4 




2.02 


0.076 


10 6 


6.65 xl0~ 4 


6.55 xlO- 4 


6.2 xlO~ 4 



(1) According to Eqs. (16) and (17). 



leading order (LO) in the coupling constant and including the next-to- 
leading order (NLO) distribution effects as an overall factor. 

Figure 4 presents the curves for pQCD calculations of the charm produc- 
tion Z-moments carried out supposing a primary spectrum either with, or 
without, the knee. The values of Z-moments for D-mesons were extracted 
directly from Ref.[13]. The Z-moments for A+ and D s are derived by taking 
Z(A C ) ~ 0.3 Z(D) and Z(D S ) ~ 0.2 Z(D), respectively, based on values 
assumed for the corresponding cross section ratios. 

More recently Gelmini, Gondolo and Varieschi (GGV)[24] updated the 
calculation to include the full contribution of NLO predictions to the lepton 
fluxes. While TIG scales the LO cross sections by a constant factor of K = 2 
to obtain the NLO contribution, GGV evaluates explicitly the NLO com- 
ponent, as tailored by Mangano, Nason and Ridolfi[25]. At the end, the net 
calculation corresponds as to multiply the LO term by an energy dependent 
factor K. In the 10 2 to 10 11 GeV energy interval, it starts at the lowest 
energies with K — 3, decreases to around 2 for most of the intermediate 
energies, increasing slightly at the high energy extreme. However, the main 
difference between the two calculations comes from the extrapolation of the 
gluon parton distribution function, which produces higher charm cross sec- 
tions even at LO. We do not duplicate here the particular effects implied 
over the Z-moments, leaving to consider the overall changes, resulting from 
GGV approach, at the evaluation of the lepton fluxes. Despite these differ- 
ences, the atmospheric cascading routines are the same in both TIG and 
GGV calculations. A discussion of an alternative evaluation of the charm 
production based on pQCD, compared to TIG's calculation, is presented by 
Pasquali et al. [14]. 

In Section 4 we compare the resulting prompt lepton fluxes, evaluated using 
each of these models. It is interesting to include in this comparison the results 
obtained by E. Zas et al. [10]. They calculate extreme cases of charm produc- 
tion, at both low and high production rate limits. As for the high end, they 
assume a charm production cross section which is 10% of the total inelastic 
cross section (called Model-A), behaving as log 2 (s) at high energies, y/s be- 
ing the center-of-mass system (c.m.s.) energy. At the low end lies a pQCD 
model at NLO, with structure functions given by Kwiecinski-Martin-Roberts- 
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Fig. 4. Charm Z-moment components for pQCD, considering the primary spectrum 
either with or without a knee, as a function of energy. 

Stirling[26], adopting one choice of parameters that leads to relatively hard 
parton distribution functions (called Model-E). 



4 Results 



4-1 High-energy atmospheric muon flux 



A reproduction of calculated fluxes of high-energy muons reported by different 
authors is presented in Figure 5. Comparing the conventional flux, from the 
decay of pions and kaons, obtained by Lipari and TIG, with the prompt flux, 
calculated by Volkova et ai, Bugaev et ai, TIG, GGV and Zas et ai, we note 
that the cross-over of the conventional and prompt components may occur at 
energies between 2 x 10 3 GeV and 2 x 10 6 GeV. The prompt flux intensity, 
itself, spans up to four orders of magnitude! 

We remark that the lowest curves, labeled TIG, GGV and Zas (E), are cal- 
culated assuming the pQCD charm production model. Although based on the 
same framework, these calculations of the charm cross section are subjected 
to theoretical uncertainties, which arises from the possible range of charm 
quark masses, as well as the factorization and renormalization scale depen- 
dence. Moreover, different assumptions are made for the parton distribution 
functions, needed at very small parton momentum fractions, not measured at 
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accelerators. The curve labeled Bugaev is calculated within the RQPM. As an 
intrinsic-charm model, it predicts relatively hard inclusive spectra. In spite of 
that, the total inclusive cross section can be made rather large, since it de- 
pends strongly on the assumptions about the charm structure function of the 
projectile hadron. The curve labeled Volkova is obtained from a parametriza- 
tion given in Ref. [7], because the original work[9] quotes only flux ratios. The 
QGSM assumed in this calculation is considered to represent quite well open- 
charm production, and it is known to describe a wide variety of data on 
hadronic interactions. However, the model predicts a preferential production 
of certain secondary particle species which is not supported by experiment. 
In addition, predictions based on this model seem to exceed the experimental 
observations of horizontal air showers, measured by AKENO[27]. Finally, the 
calculation labeled Zas (A) represents an extreme and crude upper bound for 
the prompt fluxes, based on the assumption that a charm is produced 10% 
of the time in these high-energy interactions. Actually, it is ruled out by the 
bounds set to the charm cross section from the above mentioned analysis of 
AKENO data. 

One must be careful when analyzing differences in the curves of Figure 5. 
We call attention to the fact that not only different charm models are be- 
ing compared. The reported prompt fluxes are calculated adopting different 
atmospheric cascading routines. In order to better evaluate the observed dis- 
crepancies, we proceed the analysis by taking into account the separate effect 
of each single ingredient in the shower process. Later on, we analyze the effect 
of multiple ingredients, to evaluate how the choice of a different combination 
may affect the final result. 



4-2 Single ingredient effects 

If we calculate the fluxes of prompt neutrinos [y^ and u e ), we obtain essentially 
the same values as the prompt muon flux (see for example Refs. [7,13,24]). The 
reason is that both the parent (D or A c ) and the daughter (K or A G ) particles 
are massive compared to the leptons and the decay kinematics become blind 
to lepton family number or flavor. 

The prompt lepton flux is also essentially independent of x (for depths greater 
than a few interaction lengths), due to the fact that the main contribution to 
the high-energy lepton flux must come from the first interactions of primary 
nucleons with air nuclei, while they are still energetic enough. In addition, 
for a fixed detection level, the flux is insensitive to the zenith angle 9, up to 
the charm critical energy. This is also a consequence of the fact that the bulk 
of energetic leptons are produced high in the atmosphere and will reach the 
detector, regardless the amount of atmosphere traversed (less depth in the 
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Fig. 5. Comparison of calculated differential vertical muon fluxes at sea level, as 
reported by several authors. The fluxes are multiplied by E^. 



vertical direction, larger slant depths for showers close to the horizon). For 
energies above the critical energy, the charm particle decay length becomes 
comparable to the interaction length, and we feel the effects of its angular 
dependence, given in Eq. (10). 



Since the prompt lepton flux is almost independent of lepton flavor, detection 
depth and zenith angle, hereafter we will perform our comparisons using the 
muon-neutrino vertical flux [9 = 0°) at sea level (x ~ 1000 g/cm 2 ). 

First we investigate the effect of the primary spectrum at the top of the at- 
mosphere (Section 3.1). To do so we calculate the prompt lepton flux, keeping 
all ingredients fixed, but the primary flux. Just for comparison purposes, we 
consider ingredients for the showering process used by Bugaev. The effect of 
the primary spectrum, alone, to the final result is considerable, as shown in 
Figure 6. The spread on the resulting fluxes generally increases with energy. 
For example, at 10 9 GeV, the difference between using Lipari's single slope 
and TIG's primary is a factor 10 times, while they started together at 10 3 
GeV. Also a big shift is present for the curve with AKENO primary against 
all others with knee, above 10 5 GeV. 

The next two ingredients considered are the nucleonic interaction length and 
the nucleonic Z-moment (Section 3.2). Apart from assuming a constant inter- 
action length (as done by Lipari), for which the overestimated value at high 
energies pulls the neutrino flux down, the resulting fluxes are rather insensitive 
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Fig. 6. Comparison of prompt neutrino fluxes for different primary spectrum models 
(see Table 2), assuming all the other ingredients fixed. 

to the choice of A at. Similar situation occurs with the fluxes calculated chang- 
ing only the nucleonic Z-moments, the difference being that the overestimated 
flux comes from TIG (Znn(i) ~ 0-5, at energies below 10 5 GeV), pushing the 
neutrino flux up, while other models (Z N n(i) ~ 0.2 — 0.3), result basically in 
the same final prompt flux. 

The fluxes are also rather insensitive to the charm interaction length up to 
10 7 GeV, as they should, since that is about the value of e cr iu c for charmed 
particles (see Table 1). Above this range we discriminate the models up to a 
factor of two times, whether we use Aj constant or with a log(E) dependence. 

When evaluating the effects of the choice of charm production model (Sec- 
tion 3.3), the big uncertainties in the inclusive cross sections of charm pro- 
duction are transmitted to the calculated prompt lepton fluxes, as seen in 
Figure 7. The spread between the prompt flux calculated with RQPM-SV and 
pQCD reaches a multiplicative factor of 20 at higher energies, solely due to 
the choice of Z Ni (j). 



4-3 Extreme ingredient combinations 

In the previous section the effect of each individual ingredient of the atmo- 
spheric showering process has been appreciated. Now, we investigate how the 
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Fig. 7. Comparison of prompt neutrino fluxes for different charm Z-Moment models 
(see Section 3.3), assuming all the other ingredients fixed. 

choice of multiple ingredients combined in different ways affects the resulting 
flux. Let's consider an example. The choice of the production model QGSM 
makes the flux dominates over other choices of charm production only up to 
about 2 x 10 5 GeV, if all other parameters are fixed (Figure 7). The combi- 
nation of the QGSM with the AKENO primary flux, which solely contributes 
to larger prompt flux up to 10 7 GeV (Figure 6), makes this new set dom- 
inate upon a broader range of energy. Choosing a charm interaction length 
model that is responsible for increasing the prompt flux at higher energies, we 
build a set which dominates over yet a larger range. Such a combination of 
showering parameters leads to a maximum flux configuration. With a similar 
procedure we can build a minimum configuration. In this case it is interest- 
ing to consider the effect of Lipari's primary up to the energy of the knee, 
above which Bugaev's Model-D for the primary spectrum produces the min- 
imum prompt flux. The other showering parameters may also be chosen as 
to stress the maximum or minimum possible behavior of a given flux calcu- 
lation. Table 4 presents our suggestion of multiple ingredient combinations 
selected to obtain extreme outputs. The prompt lepton flux curve calculated 
for each charm production model (QGSM, RQPM or pQCD), can be shifted 
up or down, depending on the chosen combination of ingredients, an effect 
illustrated in Figures 8 to 10. The band in each figure reflects the freedom to 
change the resulting flux between the maximum (MAX) and minimum (MIN) 
extreme combinations listed in Table 4. 
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Fig. 8. Prompt neutrino flux calculated with QGSM. The band represents the range 
between extreme ingredient combinations. Prompt flux from Volkova et al. and 
conventional neutrino fluxes also shown for illustration. 
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Fig. 9. Prompt neutrino flux calculated with RQPM. The band represents the range 
between extreme ingredient combinations. Prompt flux from Bugaev et al. and con- 
ventional neutrino fluxes also shown for illustration. 
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Table 4 

Extreme ingredient combinations of atmospheric showering parameters, selected to 
obtain maximum or minimum behavior of prompt flux. 



Ingredient 



MAXIMUM Set 



MINIMUM Set 



iV ,7 (for E < E kne 
iVo,7 (for E>E k 
X N (E) 

Znn{i) 
HE) 



AKENO Lipari 
AKENO Bugaev-D 
power-law[16] constant [12] 

scaling violation[13] constant with knee[18] 
constant [20] \og{E)[7] 
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Fig. 10. Prompt neutrino flux calculated with pQCD. The band represents the range 
between extreme ingredient combinations. Prompt flux from TIG and conventional 
neutrino fluxes also shown for illustration. 



5 Conclusion 



The calculation of the prompt lepton flux produced in the atmosphere by the 
semileptonic decay of charmed particles is rather straightforward, but we can- 
not say the same for the analysis of the results. For instance, the lack of precise 
information on high-energy charm inclusive cross-section in hadron-nucleus 
collisions is accompanied by a variety of options for the particle showering 
process in the atmosphere. We described different ingredients of the calcula- 
tion, comparing side by side several parametrizations for each one of them, and 
evaluated their relative importance to the final result. The major effects are 
due to the choice of the primary spectrum at the top of the atmosphere and, 
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of course, to the choice of the charm particle production model. Only first nu- 
cleonic interactions play essential role in determining the prompt lepton flux 
down at sea level, therefore variations in the nucleonic attenuation lengths 
are not that relevant, while the charm interaction length have some influence, 
above the charm particle critical energy. We observed how different combina- 
tions of ingredients can shift the resulting flux curves up or down. Therefore, 
the comparison between different calculations must be carried out with great 
care. The bands displayed in Figures 8-10 correspond to the allowed regions 
for the prompt lepton flux calculated respectively with the charm production 
models QGSM, RQPM and pQCD, if the showering process parameters are 
mixed differently. 

The prompt lepton crossover energy, that is the energy above which the charm 
particle decay products dominate over the conventional pion and kaon decay 
induced fluxes, is yet an uncertain quantity. According to Figures 8-10, it 
may be anywhere between 10 4 and 10 6 GeV. If the cross-over is "low enough" 
(about 10 4 GeV), then neutrino telescopes now operational can therefore take 
advantage of the isotropy of the prompt lepton flux, to search for an zenith 
angle independent component in their data. This can also be pursued by the 
analysis of the more copiously detected down-going muons. 

Exploiting the case of tau-neutrinos, which may produce a clear signature 
in high-energy neutrino detectors [28], will be addressed in a future analysis. 
There is also the need for a more comprehensive description of the available 
data on charm production cross section, and its extrapolation to higher ener- 
gies. 

The prompt lepton flux is on the order-of-the-day of operating high-energy 
neutrino telescopes, because of the background it represents. Proposed exper- 
iments, like IceCube [29], may turn the arguments the other way around, for 
their measurements with enhanced sensitivity may provide outstanding infor- 
mation on heavy quark interactions, just by discriminating atmospheric from 
cosmic neutrinos, at energies above tens of TeV. 
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